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■ Abstract. All numerical simulations of compressible convection reveal 

the presence of strong downwards directed flows. Thanks to helioseis- 
mology, such plumes have now been detected also at the top of the solar 
convection zone, on supergranular scales. Their properties may be crudely 
described by adopting Taylor's turbulent entrainment hypothesis, whose 
validity is well established under various conditions. Using this model, one 
finds that the strong density stratification does not prevent the plumes 

. from traversing the whole convection zone, and that they carry upwards 

^> ' a net energy flux. They penetrate to some extent in the adjacent sta- 

. ble region, where they establish an almost adiabatic stratification when 

I there is little radiative diffusion. These plumes have a strong impact on 

Q>^ . the dynamics of stellar convection zones, and they play probably a key 

' role in the dynamo mechanism. 

o\ : 
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a • 1. Introduction 

> '■ 

•'— j . One of the main weaknesses of stellar physics remains our poor description of 

\ thermal convection. It is true that the widely used mixing-length treatment 

■ permits to construct models which represent fairly well the gross properties of 

stars, but it fails when one attempts to apply it to more subtle processes, such 
as convective penetration, the differential rotation of the Sun, or its magnetic 
activity. 

The situation is rapidly changing however. Over the two past decades signif- 
icant progress has been achieved through numerical simulations of increasingly 
'turbulent' convection in a stratified medium. These have shown that fully com- 
pressible convection is highly intermittent, displaying strong, long-lived, down- 
wards directed flows, which contrast with the slower, random upward motions 
(Hurlburt et al. 1986; Cattaneo et al. 1991; Nordlund et al. 1992; Muthsam 
et al. 1995; Singh et al. 1995; Brummel et al. 1996; Stein & Nordlund 1998). 
These coherent structures are called plumes ('panaches' in french), by analogy 
with those observed in the Earth atmosphere. They originate in the upper 
boundary layer, where they are initiated by the strong temperature and density 
fluctuations, which arise there in the steep superadiabatic gradient. 



1 



^ -Vi, — ■'JJ '-IT 












i ^, ^ -in . 








- ' '1- - \ 
'" . I V - . J ■ . - 











Figure 1. Time-distance tomography with the MDI instrument on board 
of SOHO: the field shown measures 8 Mm in depth and 15 Mm in width. 
( http://sohowww.nascom.nasa.gov/gallery/MDI/mdi009.gi] ) 

True, these numerical results are still difficult to apply as such to the con- 
vection zone of a star, because of the huge gap in the relevant control parameters, 
the Reynolds and the Prandtl numbers. But they suggest that coherent struc- 
tures may play an important role in the dynamics of a turbulent convective layer, 
and thus open the possibility of a different description of stellar convection. 

And indeed these plumes have now been detected in the upper layers of 
the solar convection zone, thanks to time-distance tomography. This powerful 
helioseismic technique reveals the presence of downdrafts apparently associated 
with the supergranular network, in which the temperature is strongly correlated 
with the vertical velocity (D'Silva et al. 1996; Duvall et al. 1997). Fig. || shows 
an example of such tomographies; it can been downloaded from the SOHO 
gallery on the Web. 

Turbulent plumes are commonly observed in the Earth atmosphere, where 
they rise above concentrated heat sources (chimneys, nuclear plants, volcanos, 
etc.). Their theoretical interpretation is still based on the entrainment hypoth- 
esis, which was first proposed by G.I. Taylor (cf. Morton et al. 1956); we recall 
it below. In astrophysics, plumes have been invoked by Schmitt et al. (1984) as 
responsible for the penetration into the stable domain beneath the solar convec- 
tion zone. Somewhat later an estimate of that penetration has been given by 
Zahn (1991), which is also based on a crude plume model. 

More recently, Rieutord and Zahn (1995) moved a step forwards, by letting 
the plumes develop through the whole convection zone. They adapted to a 
medium which is highly stratified, in density and temperature, the treatment 
which is used in geophysical fluid dynamics. Furthermore, they took into account 
the backflow which is induced by the confinement of the plumes in a limited 
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Figure 2. The plume flow in an adiabatic atmosphere (after Rieutord & 
Zahn 1995; courtesy A&A). 

volume, and which determines the depth the plumes may reach. This review 
recalls the main properties of that model; it will also address the role the plumes 
may play in convective penetration and in the generation of gravity waves. 

2. One single plume in an isentropic stratification 

We shall first examine the basic properties of a turbulent plume which resides in 
an isentropic, plane-parallel atmosphere, assuming that the plume is fed steadily 
by a cold layer of fluid lying on the top of the atmosphere. A schematic view of 
the flow in such a plume is given in flg. ??. 

2.1. Governing equations 

The equations governing the structure of turbulent plumes have been established 
in the original work of Morton et al. (1956). One assumes that the flow is 
stationary, that the plumes are axisymmetric about the vertical, and that all 
horizontal variations, namely those of the vertical velocity Vz, of the density 
excess 6p and of the excess of specific enthalpy dh, have the same gaussian 
profile: 

Vz{r,z) = F(z) cxp(— ,s^/5^) 

6p{r, z) = Ap{z) cxp(-sV^^) (1) 
6hlr,z) = Ah{z) eicp{- s^/b'^) 

where b = b{z) is the effective radius of the plume. For convenience, we take the 
vertical coordinate z pointing downward; the density and the specific enthalpy 
of the isentropic atmosphere will be designated by po{z) and ho{z). 

The equations describing the vertical profile of the plume arc derived from 
the three basic equations of fluid mechanics expressing the conservation of mass, 
momentum and energy. 
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Conservation of mass. In the steady flow of the plume, the mass conservation 
imphes 

div{pv) = {spvs) + ^ (pvz) = ; 
when integrated over s, this equation becomes 

pv.sds + [spvs]^'^ = 0. (2) 



d_ 
dz 



The entrainment hypothesis made by G.I. Taylor postulates that the radial 
inflow of matter into the plume is proportional to the central vertical velocity: 



lim svg 

s—*+oo 



-ah{z)\vz{0, z) 



(3) 



where a is the entrainment constant. The absolute value guarantees that the 
plume is always accreting matter, whether it is directed upwards or downwards. 
Using the gaussian profile of Vz and 5p, one casts the mass equation in its final 
form 



d_ 

dz 



pob^V =2apobV, 



(4) 



in the limit of vanishing density contrast Ap/pQ <C 1. 



Conservation of momentum. Ignoring the viscous stresses, the flow obeys the 
steady Euler equation: 

djipvjVi + SpSij) = 5pgi , 

where the static equilibrium has been subtracted. We assume that the plume is 
in pressure balance with the surrounding medium; it leads us to 

\^^{spvsvz) + ^^{pvt)=5pg 

{g is the gravity, which we assume constant). After integration as before, we get 

d 



dz 



pob'V^ =2gb'Ap. 



(5) 



Conservation of energy. We start from the steady energy equation 

div(pe'i;) = div(xVr) — pdivv , 

where e is the specific internal energy and x the radiative conductivity, and 
rewrite it by using the momentum equation and the definition of the specific 
enthalpy dh = de + d{p/p): 



div 



{h+^v'')pv-xVT 



pvg. 



This equation may be further simplified by remembering that 

g = V/io 
4 



in an isentropic atmosphere (see below) . Finally, the conservation of the energy- 
flux is expressed in conservative form: 




(6) 



where we recognize respectively the enthalpy, kinetic and radiative fluxes. 

To make contact with previous work, we first neglect the radiative flux. 
Proceeding as before, we integrate in the radial direction, and obtain 



where F designates the total flux (enthalpy plus kinetic) carried by one plume. 

Prom now on, we shall drop the subscript q, since there will be no ambi- 
guity anymore between the local values of density and temperature, and their 
counterparts in the isentropic atmosphere. 

Thermodynamics. Por simplicity, we shall assume that the fluid is a perfect 
monatomic gas, with 7 = 5/3 being the adiabatic exponent. We thus complete 
the system of governing equations with the two thermodynamic relations: 



where the latter expresses the pressure balance between the plume and the 
surrounding medium (cp is the heat capacity at constant pressure). 

The isentropic atmosphere. As is well known, the isentropic atmosphere is a 
polytrope of index = 1/(7 — 1). In plane-parallel geometry, the density and 
the temperature vary with depth z as: 



The origin of the vertical coordinate z is taken at the 'surface', where pressure, 
density and temperature all vanish, and the subscript j designates the refer- 
ence level for which we choose here the base of the atmosphere. The reference 
temperature Tj is related to the reference depth by Zi = CpTi/g. 

2.2. Asymptotic behavior 

The equations governing the plume flow arc nonlinear and in most cases it is 
necessary to integrate them numerically, with boundary conditions imposed at 
the start of the plume, for instance on Ap, b and V. However, the system 
has asymptotic solutions which can easily be obtained in analytic form. The 
one which presents most interest corresponds to the developed phase of the 




(7) 



Ah = CpAT 





(8) 
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downwards directed plume, when it is fully controlled by entrainment. The 
solution then obeys the power laws 



h = h{z/z,f, Y = V,{zlzi)\ ^p = ^p,{z/z,r (9) 
where the exponents take the values: 

p = \ 

= _1±1 = 27-1 

3 " 3(7-1) (10) 

2q + 7 77-5 



m 



3(7 - 1) 



A direct consequence of Taylor's entrainment hypothesis (^) is that the plume 
radius increases linearly with depth. The half-angle (3 =aXaD.[hi / Zi) of the cone 
is related to the adiabatic exponent: 

3a 3a(7 — 1) ,^ ^ , 

tan/3= = — ^- (11) 

g + 2 27-1 ' 

With 7 = 5/3 (perfect monatomic gas), we have tan/3 = 6a/7, whereas tan/3 = 
6q/5 in the absence of stratification. The entrainment coefficient a has been 
determined experimentally in various situations, and the value 0.083 is widely 
adopted (Turner 1986); we shall assume that this value also applies to our fully 
ionized gas. 

For given initial conditions imposed at the top of the plume, the solutions 
rapidly reach their asymptotic regime: the plumes start contracting while they 
accelerate, and then they take their characteristic conical shape (cf. Rieutord & 
Zahn 1995). 

Note that one does not retrieve the incompressible case by letting the poly- 
tropic index tend to zero, g — > 0. The reason for this singular limit is the assump- 
tion made in the classical treatment (Morton et al. 1956) that the surrounding 
medium is isothermal, while in our stratified atmosphere the temperature grows 
linearly with depth. Let us emphasize however that our asymptotic solution (P) 
is rather special, since its focal point z = coincides with the 'surface' of the 
atmosphere. 

2.3. Turbulent plumes carry upwards a net energy flux 

An interesting property of this asymptotic regime is the strict proportionality 
between the flux of kinetic energy and that of enthalpy: 

f -pV^sds 

R = = est. (12) 

pep 6T Vsds 



With the gaussian profile we have adopted for the plume, this ratio is given by 

R = -^ = -'-2^. (13) 
g + 2 27-1 ^ ^ 
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Its value is 1 for the Boussinesq fluid (no density stratification) but it is 4/7 
in the isentropic atmosphere (for the perfect monatomic gas), meaning that the 
net energy flux is directed upwards. 

In their numerical experiment Cattaneo et al. (1991) observed to their sur- 
prise that the strong downdrafts transported a flux of kinetic energy which very 
nearly canceled the flux of enthalpy. The reason is that their plumes were not 
turbulent, but laminar flows which satisfied approximately Bernoulli's theorem. 
In recent simulations, which have been achieved with better spatial resolution, 
there is more entrainment, and the plumes do carry a net energy fiux upwards. 



3. Towards a more realistic model; application to the solar convec- 
tion zone 

Up to now, we have considered the case of a single plume in an unlimited isen- 
tropic plane layer. To progress towards a more realistic situation, we must take 
into account the spherical geometry of the fluid layer and include the variations 
of gravity with depth. Here we shall apply the model to the solar convection 
zone (SCZ); the convective core with its rising plumes has been treated by Lo 
and Schatzman (1997). 

We neglect the mass of the envelope, so that the density and temperature 
field of the background are given by 

Tiz) = T, (--I 



\ — T] \X 

q 



where x = r/R^ is the radial coordinate scaled by the radius, and rj = Vcz/Rq 
its value at the base of the SCZ (our reference level). 

We shall also take into account the radiative fiux in the energy equation, 
thus (0) now reads 

^pb^VAh + ^pb^V' = -(Ftot - Frad)/vr . (15) 

With our assumption that all the convective energy is transported by the plumes, 
whose number is A^, the total flux amounts to the luminosity of the star: 

NFtot = LQ. (16) 

The radiative contribution to the flux carried by one plume is Fradj where 

2 fdT\ 16aT^ 



\dz J^d Skp 
with X being the radiative conductivity, k the opacity, and a the Stefan constant. 
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3.1. Plumes diving in a rising counterflow 

We shall now take into account that there are numerous plumes and that they 
share a finite volume. This has two consequences: the plumes will interact and 
possibly merge, and they will move in a rising counterflow. Let us first estimate 
the maximum number of plumes which may coexist in the SCZ, assuming that 
they all originate at the top and that all reach its bottom without being hindered 
by the counterflow. This number is given by the ratio of the area of the base of 
the SCZ to the section of a single plume at this depth: 



N 



7] 



(1 — T]) tan P 



where /3 is given by (pTl), and i] = Tcz/Rq is the scaled radius of the base of 
the SCZ. Taking r] = 0.7, the result is A'max ~ 4300. We see that even if the 
plumes were closely packed, they would not be very numerous. But they will be 
actually less, because as soon as their occupy more than half of the total area, 
at depth 2;j/\/2, the velocity of the backflow exceeds that of the plumes, and 
this will affect their dynamics. 

To model the effect of this counterflow, we shall assume that the entrainment 
of mass into the plume is proportional to the relative velocity of the plume with 
respect to the surrounding fluid. This yields the mass conservation equation 

-^\pb^v]=2apbiV-Vu) (17) 
dz L J 

where Vu is the upward velocity of the surrounding fluid; since < 0, entrain- 
ment is enhanced. The momentum equation needs also be completed by a term 
taking into account the entrainment of momentum into the plume. We thus 
transform (0) into 



d_ r 

dz 



p¥V' =2gb^Ap + 4apbVu{V -Vu). (18) 

Equations (|17|) and (jlj) are completed by 

(47rr2 - NTTb'^)pVu + Nnb'^pV = (19) 

which expresses the global conservation of mass. The consequence of these 
additional effects is that plumes are able to reach the bottom of the convection 
zone only if they are strong enough, and not too numerous. From the numerical 
integration of these equations it turns out that only about 1000 plumes are able 
to reach the bottom of the SCZ (see Rieutord & Zahn 1995). 

3.2. Plume coalescence 

At the top of the SCZ, the number of plumes is probably of the same order as 
that of the granules, which is considerably larger than 1000. Hence a drastic 
reduction of the plumes number must occur as depth increases. The plumes 
are stopped at a level which depends on their strength: only the 1000 strongest 
are able to reach the bottom. The others will either dissolve in the (turbulent) 
backflow and disappear as such, or they will merge with a stronger companion. 
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Such merging is observed in the numerical simulations (Stein & Nordlund 
1989; Spruit et al. 1990): two neighboring plumes are advected by each other, 
since they both accrete surrounding fluid, and they are pulled to each other until 
they finally coalesce. The merging depth can be estimated by assuming the flow 
be in the asymptotic regime; the result is 

' (20) 



where d is the initial separation of the two plumes, which we assumed here to 
be identical. Using the value of fi for the monatomic ideal gas and a = 0.083, 
we get Zm ~ 9.2 d . 

Similar coalescences will then repeat at increasing depths, until the sur- 
vivors reach the base of the convection zone, in a sort of inverse cascade with a 
large scale flow building up from smaller scales. 



4. Penetration at the base of a convection zone 

When the plumes reach the bottom of the unstable domain, they still possess 
a finite velocity which enables them to penetrate some distance into the stable, 
subadiabatic region, where they establish a nearly adiabatic stratification by 
releasing their entropy when they come to rest. A first attempt to estimate 
the extent of penetration of such plumes was made by Schmitt et al. (1984). 
However, they did not include the return flow, and they had to impose both the 
velocity V and the filling factor / at the base of the unstable domain because 
they did not solve the plume equations above that level. They found empirically 
that the penetration depth varies as /^/^y^/^, and this scaling can be easily 
explained. 

4.1. An estimate of the penetration penetration depth 

The stratification at the base of the a convective envelope is sketched in fig. (P). 
(A) designates the unstable region, where the temperature gradient is main- 
tained close to adiabatic by the convective motions. One thus assumes that the 
radiative leaks are negligible compared to the advective transport of heat, which 
is the case when the Peclet number is much larger than unity: 

Pec=-»1, K = -^, (21) 
K pCp 

where K is the thermal diffusivity and x the radiative conductivity, v and £ are 
the velocity and size characterizing the convective motions, i.e. here the central 
velocity V and the width b of our plumes. At the base of the SCZ, Pec ~ 10^. 

Due to the steady increase with depth of the radiative conductivity X; the 
radiative flux F,.^ = x{dT /dz)f^ rises until it equals the total flux Ftotai at the 
level z = Zj, where also the radiative gradient equals the adiabatic gradient: 
(dr/(iz)rad = {dT/dz)ad- If there were no convective penetration, this would 
be the edge of the convection zone, as predicted by the Schwarzschild criterion, 
and the temperature gradient would thereafter decrease as dT/dz = -Ftotai/x- 
But the motions penetrate into the stable region (B) and they render it nearly 
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Figure 3. Schematic structure of the base of a convective envelope (after 
Zahn 1991; courtesy A&A). 



adiabatic over some distance dpen, while being decelerated by the buoyancy force. 
When the Peclet number has dropped below unity, the temperature gradient 
settles from adiabatic to radiative in a thermal boundary layer (C). 

Since the penetration depth is rather small, one may simplify the problem by 
ignoring the variation with depth of most quantities (density, width of the plume, 
etc.), and by neglecting the kinetic energy flux and the turbulent entrainment 
(or detrainment). One keeps of course the variation with z of the conductivity; 
in the vicinity of Zi, the radiative flux is approximated by 



fdT\ 

-trad — X -y- — -Ttotal 
V / ad 

Therefore the convective flux varies as 



{z-z, 

az 



(22) 



-^conv = --^total ( 3^ ) {z - Zi) . (23) 



dz 



This enthalpy flux may be expressed in terms of the vertical velocity and of the 
temperature contrast in the plumes, as in (0): 

F,on. = -fpCpVAT, (24) 

where we have introduced the filling factor / of the plumes, defined as the 
fractional area covered by them: / = Nb'^ /R'^^. 

To estimate the penetration depth, we follow the plumes from z = Zi, 
where their velocity is Vi, until they stop at z = Zi + dpen- Their deceleration is 
described to first order by 

IdV^ Ap AT 

0-^=9 — = -g—. (25) 
2 dz p 1 
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After elimination of AT with (p^, the integration of ( pSD yields the following 
expression for the penetration depth dpen: 



'pen 




) 



-HpH^Q 



(26) 



where Hp is the scale-height of the pressure and that of the radiative con- 
ductivity. This is the scaling obtained empirically by Schmitt et al. (1984). 

At that point one may argue that the term Q in curly brackets is determined 
by the dynamics in the convection zone, and that it will not depend much on 
its depth, provided this depth is large enough and that convection is sufficiently 
adiabatic. Q is probably smaller than unity: in the mixing length treatment, 
Q = (1/10) A/ Hp in the bulk of the convection zone. This figure, with A/ Hp ~ 
1.5 and Hy, Hp/2, would yield a penetration depth of the order of 1/5 of the 
pressure scale-height at the base of the solar convection zone. 

4.2. Penetration of plumes originating at the top of the convection 



Rieutord and Zahn (1995) computed the penetration in a more consistent way, 
by integrating the governing equations from the top of the convection zone until 
the plumes come to rest. When neglecting the backflow, the penetration extends 
to about one half of the pressure scale-height, independent on the number of 
plumes. But the picture changes drastically when one includes the interaction 
of the plumes with the backflow. Then the penetration depth depends on the 
number of plumes: it decreases almost linearly with A'^, from one pressure scale- 
height for A — > to for A ~ 1200. The reason for this trend is that the 
backflow increases with the number of plumes, and that it slows down their 
motion by loading them with upward momentum. 

Unfortunately, this model is unable to predict the number of plumes. More- 
over, a much more sophisticated model is needed to describe the termination of 
the plume flow, which would include the loss of mass by the plume in this region, 
or 'detrainment', and its transfer to the upward motion. 

4.3. Penetration vs. overshooting 

What we have called convective penetration is that part of the excursion of 
convective motions into the stable region below where they enforce an almost 
adiabatic stratification (Vad — V ^ 1). It corresponds to region (B) in fig. ^, in 
which the Peclet number is substantially larger than unity (we recall that at the 
base of the SCZ, Pec ~ 10^). In comparison, the thermal adjustment layer (C) 
is very thin: about 1 km in the Sun, assuming that all plumes have the same 
strength (Zahn 1991). 

The situation is very different in stars which possess a shallow convection 
zone. For instance, at the base of the convection zone of an A-type star. Pec ~ 
1/100, which means that the convective motions which enter the stable layer 
below rapidly settle in thermal equilibrium with their surroundings, and that 
they hardly feel the buoyancy force (Toomre et al. 1976; Freytag et al. 1996). 
The region (B) has disappeared - all what remains is region (C). To insist on 



zone 
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the contrast with the almost adiabatic penetration at high Peclet number, we 
prefer to call this convective overshoot (or undershoot). 

This is not only semantic matter. To take one example, numerical simula- 
tions which are performed to represent A-type stars cannot be used to predict 
the amount of penetration below the solar convection zone, as was attempted 
by Blocker et al. (1998). In fact, for lack of spatial resolution, present day simu- 
lations are unable to achieve a Peclet number which would realistically describe 
convective penetration; for instance, in the calculations made by Hurlburt et al. 
(1994), regions (B) and (C) are of comparable size, and a substantial fraction of 
the extension into the stable region is due to overshoot. 

4.4. Generation of internal gravity waves 

In the laboratory, convective penetration is observed to generate internal waves, 
which propagate in the stable adjacent layer (Townsend 1966; Adrian 1975). the 
same is also expected at the boundary of stellar convection zones. Such waves 
have indeed been observed in numerical simulations of penetrative convection, 
where they seem to be produced by the pronounced downdrafts (Hurlburt et 
al. 1986; Hurlburt et al. 1994), which are the 2-dimensional analogues of 3- 
dimensional plumes. In these calculations one observes a strong feedback of 
the waves on the downdrafts, because the waves are reflected back by the lower 
boundary of the computational domain, a property which is not expected from 
the solar gravity waves produced by the plumes. 

The generation of gravity waves by turbulent plumes is presently being in- 
vestigated by M. Kiraga, in 2- and 3-dimensional simulations. He takes great 
care of avoiding the reflexion of the waves by the lower boundary of his compu- 
tational domain, by implementing in its vicinity a strong viscous damping. The 
first results have been presented at the Granada workshop on convection (Ki- 
raga et al. 1999); they look promising, and will be used to test the prescriptions 
which have been proposed for estimating the flux of gravity waves. 

Our interest in these gravity waves is motivated by the fact that they may 
extract angular momentum from the radiative interior of solar-type stars (see 
Kumar & Quataert 1997; Zahn et al. 1997; Kumar et al. 1999). 

5. Summarizing the properties of turbulent plumes 

There is little doubt that the convection zone of solar-like stars exhibits strong 
spatial intermittency, with downwards directed plumes carrying a substantial 

fraction of the energy flux. Such plumes are seen in all numerical simulations, 
and they have now been observed in the uppermost part of the solar convection 
zone. 

If Taylor's parametrization of the entrainment of surrounding fluid into 
turbulent plumes holds in the hot plasma of stellar interiors, and we see no 
reason why it should not, then we expect such plumes to traverse the whole 
convection zone. According to the crude model presented above in §2, these 
plumes behave much as if they were traversing an isentropic atmosphere, with 
no feed-back at all, and ignoring the radiation flux. Shortly after their start, 
they reach an asymptotic regime with their size increasing linearly with depth. 
Due to the density stratiflcation, the cone angle /3 is somewhat smaller than in 
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the Boussinesq case treated so far: tan/3 = 6q;/7 instead of 6q;/5, with a being 
the entrainment constant. 

These plumes carry kinetic energy downwards, i.e. in the 'wrong' direction. 
But their enthalpy flux always exceeds the kinetic flux: the ratio between the 
two is constant in the asymptotic regime, and its value is 4/7, if one assumes 
that the horizontal profile of the plume is gaussian. 

The maximum number of plumes that may reach the base of the solar 
convection zone, when taking the counterflow into account, is around 1000, 
which means that the spacing between plumes is about 60 Mm. Note however 
that this evaluation is based on a very strong assumption: namely that the 
totality of the convective energy flux (enthalpy plus kinetic) is carried by the 
plumes, thus neglecting the contribution of the interstitial medium. 

Finally, under the same assumption, the extent of penetration, at the base 
of the convection zone, depends on the number of plumes N reaching that level: 
it varies from for A'^ ~ 12000 to one pressure scale-height for N ^ 0. Un- 
fortunately, the number of plumes cannot be predicted with the simple model 
presented above. Let us recall that helioscismology yields a value of about O.lHp, 
with the assumption that there is sharp transition from adiabatic to radiative 
slope at the edge of the penetration layer (Roxburgh &; Vorontsov 1994). With 
plumes of different strength, the transition would be smoother and more difficult 
to diagnostic. 

All these results are independent of the top boundary conditions on the mass 
flux and on the initial momentum flux carried by the plumes. What determines 
the flow is the convective energy flux near the surface, which would be conserved 
along the plume if there were no radiative transport. Unfortunately the simple 
model presented above is unable to render the complex dynamics of these layers, 
where important radiative exchanges occur. They can only been described by 
full-fledged 3-dimensional hydrodynamical simulations, such as performed by A. 
Nordlund, R. Stein and their collaborators. However we cannot rule out that a 
sophisticated closure model, well beyond what has been described above, may 
succeed in better capturing the main properties of such turbulent plumes (see 
V. Canuto's contribution to this workshop). 

How do these plumes change our picture of stellar convection zones? Plumes 
presumably exert an important effect on the thermal structure - that is, on the 
small deviations from the isentropic stratification which are known to induce 
meridional circulation. The reason for this is that the heat transport by the 
plumes is of advective nature, whereas in most models so far, as for example in 
the mixing-length treatment, it has been represented by a diffusive process. 

They will also play a major role in the dynamics of a convection zone, and 
in determining its rotation profile. Since the ambiant angular momentum will be 
entrained into the plumes, these will spin up with depth, and they will concen- 
trate both vorticity and helicity. Owing to their robustness, such vortices could 
well be the cause of the observed differential rotation, especially of its peculiar 
property of being quasi-independent of depth (Brown et al. 1989; Kosovichev et 
al. 1997). 

Furthermore, one can easily imagine that plumes play a major role in the 
generation and in the confinement of magnetic fields. Since they are sites of high 
helicity, they will be the seat of strong alpha-effect. The dynamo simulations of 
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Brandenburg et al. (1996) have shown that laminar diving plumes are indeed 
playing an important part in the field generation and its storage at the base of the 
convective layer. In more recent simulations, they have been clearly observed to 
'pump down' the magnetic field, which is definitely an advective process (Dorch 
1998). 

To conclude, we have yet to assess all the consequences of the presence of 
plumes in stellar convection zones. For that we have to rely mainly on numer- 
ical simulations, since the stratification is much stronger in stars than in the 
Earth atmosphere, and since it cannot be reproduced in the laboratory. These 
simulations have not reached yet a sufficient resolution to adequately describe 
turbulent plumes, but they pace, and their results can be 

checked through helioseismology. 
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